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We investigate the phase behavior of a system of hard equilateral and right-angled triangles in 
two dimensions using Monte Carlo simulations. Hard equilateral triangles undergo a continuous 
isotropic-triatic liquid crystal phase transition at packing fraction <f> = 0.7. Similarly, hard right- 
angled isosceles triangles exhibit a first-order phase transition from an isotropic fluid phase to a 
rhombic liquid crystal phase with a coexistence region f> £ [0.733,0.782]. Both these liquid crystal 
phases undergo a continuous phase transition to their respective close-packed crystal structures 
at high pressures. Although the particles and their close-packed crystals are both achiral, the 
solid phases of equilateral and right-angled triangles exhibit spontaneous chiral symmetry breaking 
at sufficiently high packing fractions. The colloidal triangles rotate either in clockwise or anti¬ 
clockwise direction with respect to one of the lattice vectors for packing fractions higher than <fr x . 

As a consequence, these triangles spontaneously form a regular lattice of left- or right-handed chiral 
holes which are surrounded by six triangles in the case of equilateral triangles and four or eight 
triangles for right-angled triangles. Moreover, our simulations show a spontaneous entropy-driven 
demixing transition of the right- and left-handed “enantiomers”. 


Chirality plays an important role in nature, chemistry, 
and materials science. An object is chiral if it is not 
identical to its mirror image. The most well-known ex¬ 
ample of a chiral object is the human hand, where the left 
hand cannot be superimposed on its mirror image, the 
right hand. Also many biologically active molecules are 
chiral, e.g., amino acids are left-handed, whereas sugars 
are right-handed. The microscopic chirality of the con¬ 
stituent particles may subsequently lead to a macroscopic 
chirality of the self-assembled higher-ordered structures, 
e.g., left-handed amino acids form right-handed helical 
protein structures, and right-handed sugars lead to right- 
handed DNA double helices. Additionally, chirality is 
present in so-called cholesteric phases, which are nematic 
liquid crystals with a helical structure of the director field 
and which are frequently used in optoelectronic appli¬ 
cations, such as liquid crystal displays of laptop com¬ 
puters, cell phones, and flat screen televisions pQ. Re¬ 
cently, chiral nanostructured materials have also received 
much attention due to their intriguing optical properties 
such as a huge optical activity, strong circular dichro- 
ism, photonic band gaps, and negative refractive indices 
W- However, despite the huge amount of work de¬ 
voted to chirality, the underlying microscopic features of 
the building blocks responsible for the formation of chi¬ 
ral self-assembled structures is extremely subtle and not 
well-understood. 

Even the most basic question if particle shape alone 
can lead to macroscopic chiral structures is still unknown. 
For instance, it has been theoretically demonstrated that 
an entropy-driven isotropic-cholesteric phase transition 
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exists for hard helical particles, but these predictions 
have never been verified experimentally or by computer 
simulations [SHE], a more intriguing question would be 
whether or not achiral particles can self-assemble into 
chiral structures. Very recent experiments by Mason et 
al. on equilateral triangular colloidal platelets show an 
entropy-driven phase transition from the isotropic liquid 
to a triatic liquid crystal phase that displays three-fold 
symmetric orientational order j9]. Surprisingly, at suffi¬ 
ciently high densities, small domains of chiral dimer pairs 
that are laterally shifted in one or the opposite direc¬ 
tion, appear spontaneously in the triatic phase. The au¬ 
thors conjectured that the spontaneous local chiral sym¬ 
metry breaking is due to an increase in rotational entropy 
and may be explained by a simple rotational cage model 
El HD]. However, a recent simulation study explained the 
emergent chirality observed in these experiments by the 
rounded corners of the particles which lead trivially to 
two degenerate crystal lattices of chiral dimer pairs at 
close-packing, thereby casting doubts on the role of rota¬ 
tional entropy on the chiral symmetry breaking m- In 
addition, these simulations showed that the chiral sym¬ 
metry breaking is absent for perfect triangles, i.e ., no 
particle corner rounding, which is to be expected as the 
close-packed structure of perfect triangles is an achiral 
triangular lattice. These findings are also consistent with 
a previous simulation study on perfect equilateral tri¬ 
angles, which shows only a simple transition from the 
isotropic to a liquid crystal phase at packing fraction 
(j) = Na p /A = 0.57 with N the number of particles, A 
the area of the simulation box, and a p the particle area 

m- 

In this paper, we reexamine the phase behavior of 
hard equilateral triangles by extensive Monte Carlo sim- 
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Figure 1. Candidate close-packed crystal structures: (a) 
Equilateral triangles with two particles in the unit cell forming 
a hexagonal dimer lattice or a triatic crystal. Right-angled 
triangles with a rhombic lattice with two and four particles in 
the unit cell in (b) and (c), respectively. We show four unit 
cells for all the candidate close-packed crystal structures and 
we used red to indicate a single unit cell. 


ulations. Surprisingly, we find the spontaneous forma¬ 
tion of a novel chiral crystal phase, where the individ¬ 
ual particles spontaneously undergo either a clockwise 
or anti-clockwise rotation with respect to one of the lat¬ 
tice vectors which give rise to a regular lattice of anti¬ 
clockwise or clockwise chiral holes. We find a similar 
chiral crystal phase in a system of right-angled trian¬ 
gles. More surprisingly, we also observe a spontaneous 
entropy-driven demixing transition of the “enantiomers” 
into left-handed and right-handed chiral phases. 


I. RESULTS 

Hard equilateral and right-angled isosceles triangles 
tile the space in infinitely many ways as the rows and the 
columns of these triangles can be shifted without affect¬ 
ing their maximum packing density. At finite pressures 
hard triangles may form liquid crystal phases with ori¬ 
entational (quasi)-long-range order or solid phases with 
orientational and translational (quasi)-long-range order. 
To determine the phase behavior of hard equilateral and 
right-angled triangles, we perform Monte Carlo (MC) 
simulations and free-energy calculations. We use the 
separating axis theorem to detect particle overlaps m- 
We perform variable-rectangular-box isothermal-isobaric 
Monte Carlo simulations O ES], in which we fix the 
number of particles, A = 3000 — 13000, the pressure P, 
and the temperature T. We compress the system from 
a low-density isotropic fluid phase to a solid phase by 
slowly increasing the pressure. We observe that a system 



Figure 2. f ex + log N/N of the two candidate crystal struc¬ 
tures for right-angled triangles as a function of 1/A at packing 
fraction </> — 0.91. Here f ex — F ex /NkBT = (F — F^/NksT 
is the excess free energy per particle, F is the Helmholtz free 
energy and Fid the free energy of an ideal gas at the same 
packing fraction. We observe that the rhombic lattice with 
four particles in the unit cell has a lower free energy compared 
to the rhombic lattice with two particles in the unit cell. 


of equilateral triangles undergoes a transition from an 
isotropic fluid phase to a triangular lattice with two par¬ 
ticles in the unit cell as shown in Fig. [If- On the other 
hand, right-angled triangles never crystallized within the 
simulation times that we considered, but only small 
rhombic crystalline domains with either two particles or 
four particles in the unit cell as shown in Figs. EM 
appeared spontaneously in the system. 

In order to determine the most stable thermodynamic 
phase of the two candidate crystal structures for right- 
angled triangles, we employ the Frenkel-Ladd method as 
described in Refs. [16] jlT] to compute the free energies 
of both rhombic phases at packing fraction </> = 0.91. 
For more details regarding the implementation of this 
method, we refer the reader to Refs. ESI EH]. We show 
f ex + log N/N as a function of 1/A in Fig. [2] for both 
candidate crystal structures. Here f ex = F ex /NkBT is 
the excess free energy per particle, ks denotes Boltz¬ 
mann’s constant, A the number of particles, and T the 
temperature. We find in agreement with Ref. [T9] that 
f ex + log N/N is a linear function of 1/A with the inter¬ 
cept at 1/A = 0 corresponding to the excess free energy 
for infinite system size. If we extrapolate the excess free 
energy to the thermodynamic limit (A —f oc), we ob¬ 
serve that the rhombic lattice with four particles in the 
unit cell has a lower free energy than the one with two 
particles in the unit cell. 

Subsequently, we determine the equations of state 
(EOS) from compression runs using the isotropic fluid 
phase as initial configuration in NPT Monte Carlo sim¬ 
ulations with a variable box shape. Similarly, we obtain 
the EOS by expanding the stable close-packed crystal 
structures in NPT Monte Carlo simulations. To char¬ 
acterize the phases at high density, we determine the 
positional and orientational order at different packing 
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Figure 3. (a,b): Equations of state for hard equilateral and right-angled triangles, respectively. Both compression and 

expansion runs are obtained for a system size of N — 3200 particles for equilateral triangles and N = 1600 particles for 
right-angled triangles using NPT simulations with a rectangular box. (c): Six-fold bond-orientational (^6°) and molecular 
orientational order parameters as a function of packing fraction cp for a system of hard equilateral triangles. Both the 

order parameters show a transition around cp ~ 0.7 indicating a phase transition between the liquid and triangular crystal 
phase, (d): Eight-fold bond-orientational 2 pf° and molecular orientational ^s 10 order parameters as a function of packing 
fraction <p for right-angled triangles. The coexisting densities calculated using free energies for the right-angled triangles are 
< i> = 0.733 and 0.782, and are indicated by the dotted vertical lines. Figures (e,f) show the phase diagram for the two particle 
shapes using different colors as indicated. TLCP and RLCP represent the triatic and rhombic liquid crystal phase while T x 
and R x represent their chiral triangular and rhombic crystal structures, respectively. The white region between the fluid and 
the RLCP in Fig. §o indicates the coexistence region. 


fractions <p>. To this end, we measure the n-fold bond- 
orientational and molecular orientational order parame¬ 
ters. The n-fold bond-orientational order parameter is 
given by 


rameter reads 
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where 6ij is the angle between the vector, connecting 
particle i and its nearest neighbor j, and an arbitrary 
reference axis. Here nn = 3, is the number of nearest 
neighbors. The n-fold molecular orientational order pa- 


where 0i is the angle between particle i and a fixed ref¬ 
erence axis. Here we use the x-axis as the reference 
axis. Depending on the local symmetry of neighboring 
particles around a single particle in their corresponding 
close-packed structures we set n = 6 for equilateral tri¬ 
angles and n = 8 for right-angled triangles. We calculate 
these order parameters at varying packing fractions using 
Monte Carlo simulations of N = 12800 triangles in the 
canonical ensemble, z.e., the area A of the simulation box 
is kept fixed. Additionally, we measure the spatial corre- 
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Figure 4. The chemical potential fi/k^T as a function of the 
reduced pressure Pa p /kBT for both the isotropic fluid and 
triatic (liquid) crystal phase of equilateral triangles (a) and for 
the isotropic fluid and rhombic (liquid) crystal phase of right- 
angled triangles (b). The fluid and liquid-crystal branches 
cross for right-angled triangles indicating a first-order phase 
transition, whereas there is no crossover within the numerical 
precision of our data in the case of equilateral triangles. 


lation functions for the translational, bond-orientational 
and molecular orientational order, z.e., g(r), g§ °(r), and 
^°(r), respectively, for various packing fractions in or¬ 
der to determine whether the isotropic phase transforms 
into a liquid crystal or a crystal phase. 


A. Equilateral triangles 

We first discuss our results for equilateral triangles. In 
Figs. [§a, c), we show the equation of state (EOS) along 
with the bond orientational and molecular orientational 
order parameters as a function of packing fraction (j). 
Fig. I a) displays the EOS as obtained from both the 
compression and expansion runs. We observe that the 
system undergoes a continuous phase transition from an 
isotropic fluid phase to an ordered phase with three-fold 
symmetric orientational order upon compression. In ad¬ 


dition, we observe that the close-packed triangular crys¬ 
tal melts continuously in an isotropic fluid phase during 
our expansion runs. In Fig. j3^c) , we plot the 6—fold 
bond-orientational order parameters and molecu¬ 

lar orientational order parameters 'ipQ 10 as a function of 
packing fraction 0. Fig. |3jc) clearly shows that the sys¬ 
tems develop bond-orientational and molecular orienta¬ 
tional order for 0 > 0.7 indicating a continuous phase 
transition from an isotropic fluid to a triatic phase. We 
note that the bond order parameter value is always 
lower than that for the molecular order at all pack¬ 
ing fractions. 

In order to characterize the triatic phase in more detail, 
we also measure the correlation functions for the trans¬ 
lational, bond-orientational and molecular orientational 
order for various packing fractions around the phase tran¬ 
sition using Monte Carlo simulations of N = 12800 par¬ 
ticles in the canonical ensemble. The results are shown 
in Figs. |5ja,b,c). The radial distribution function g(r) 
which indicates the correlations in the translational order 
show exponential decay for packing fractions </> < 0.87, 
and quasi-long-range decay at 0 > 0.87, which is to 
be expected as truly long-range positional order is not 
possible in two dimensions j20H22| . The 6-fold bond- 
orientational gf°(r) and 6-fold molecular orientational 
^ 0 (r) correlation functions show quasi-long-range ori¬ 
entational order for > 0.7 within the system sizes that 
we used. The presence of (quasi) long-range bond or¬ 
der and molecular orientational order and the absence of 
long-range positional order for 0 > 0.7 are characteristic 
of liquid crystalline phases [91. Hence, we find that a sys¬ 
tem of equilateral triangles undergoes a continuous phase 
transition from an isotropic fluid phase to a triatic liquid 
crystal phase at packing fraction (j) = 0.7. Upon further 
compression, the triatic liquid crystal phase transforms 
continuously into a crystal phase at a packing fraction 
(j) > 0.87. 

To corroborate our findings, we also compute the free 
energies for equilateral triangles using the Frenkel-Ladd 
method m- We use the Widom particle insertion tech¬ 
nique to determine the chemical potential and hence the 
free energy of the isotropic fluid phase at fixed den¬ 
sity. Using thermodynamic integration of the equa¬ 
tion of states we compute the free energy per particle 
/ = F/{Nk&T) as a function of packing fraction for the 
isotropic fluid, triatic liquid crystal and triatic crystal 
phases. Subsequently, we determine the phase behav¬ 
ior. To this end, we first compute the chemical potential 
ji/ksT of both systems from the free energies and plot 
it as a function of reduced pressure Pa p /fceT as shown 
in Fig. |4ja). The fluid and liquid crystal branch do not 
cross in the case of equilateral triangles, which supports 
our finding that the isotropic fluid-triatic liquid crystal 
phase transition is continuous. In addition, we find that 
the liquid crystal branch transforms continuously into 
the solid branch, indicating a continuous triatic liquid 
crystal-triatic crystal transition. 
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Figure 5. Positional g(r), bond-orientational g^° (r) and molecular orientational g^° (r) correlation functions at varying 
packing fractions 0 as labeled for equilateral and right-angled triangles. The left column contains the correlation functions 
for equilateral triangles and the right column is for right-angled isosceles triangles. All the plots are on log-log scale. (a,d): 
radial distribution function | g(r) — 1| decays algebraically for all the packing fractions. (b,e): n —fold bond-orientational order 
correlation functions g^° (r) where n — 6 and n = 8 for equilateral and right-angled triangles respectively. (c,f) n— fold 
molecular orientational order correlation functions g^° (r) with the same values of n as above. 


B. Right-angled isosceles triangles 

We now turn our attention to the right-angled isosceles 
triangles. In Fig. 3^b), we present the equation of state 
(EOS) as obtained from both the compression and expan¬ 
sion runs. Upon compression of the isotropic fluid phase, 
we observe no crystallization during out NPT simula¬ 
tions, but only the spontaneous formation of small crys¬ 
talline domains. In addition, we observe that the rhombic 
crystal phase with four particles in the unit cell, which 
is the stable crystal phase according to our free-energy 
calculations, undergoes a first-order phase transition to 
an isotropic fluid phase at sufficiently low pressures. The 
8-fold bond orientational and molecular order parame¬ 
ters, 'ipf 0 and V’s 10 ; as displayed in Fig. show that 

the system develops bond-orientational and molecular- 
orientational order for 0 > 0.7. We note again that the 
bond order parameter value * s always lower than 


that for the molecular order 'ip^ LO for all values of </>. 

In order to investigate the range of the positional 
and orientational order of the rhombic phase, we calcu¬ 
late the correlation functions for the translational, bond- 
orientational and molecular orientational order as a func¬ 
tion of packing fraction <p. We present the correlation 
functions in Figs. [5jd,e,f). Again, we find that the g(r) 
shows exponential decay for 0 < 0.89, and becomes only 
quasi-long-range for 0 > 0.89. The bond-orientational 
gf°( r ) and molecular orientational ^ , 8 d0 (r) correlation 
functions show quasi-long-range orientational order for 
(j) > 0.79 for the system sizes that we used. We thus find 
that a system of right-angled triangles undergoes a first- 
order phase transition from an isotropic fluid phase to a 
rhombic liquid crystal phase, and shows subsequently a 
continuous phase transition to a rhombic solid phase at 
</> = 0.89. 

To determine the phase boundaries of the isotropic 
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Figure 6. Chiral symmetry breaking in the solid phase, (a): Triangular lattice with triangles that display no orientational 
displacement, resulting in an achiral crystal phase, and with triangles that are shifted clockwise, yielding a triatic solid phase 
with chiral holes. (b,c): Sign notation for anti-clockwise and clockwise orientational displacements 6 of the triangles with respect 
to a fixed lattice vector. The orientation of the triangles are denoted by an arrow. The triangles that exhibit no rotational shift 
are colored green. The particles that have an anti-clockwise orientational displacement are colored blue (+), and the particles 
with a clockwise orientational displacement are colored red (—). (d): Probability distribution of the orientational displacement 
0 of equilateral triangles at varying packing fractions 0 > 0.85 as labeled. For > 0.89, we find that P{6) shows three distinct 
peaks. 
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fluid-rhombic liquid crystal phase transition, we deter¬ 
mine the free energies of the two phases using the meth¬ 
ods as described above in Sec. IA Fig. [4jb) shows the 
chemical potential li/ksT as a function of reduced pres¬ 
sure Pap/ksT for the isotropic fluid, rhombic liquid crys¬ 
tal, and rhombic crystal phase. We find a clear crossover 
of the fluid and rhombic liquid crystal branch correspond¬ 
ing to a first-order phase transition with a coexistence re¬ 
gion G [0.733,0.782]. Additionally, the rhombic liquid 
crystal branch transforms continuously into a rhombic 
crystal branch, and hence the transition from a rhombic 


liquid crystal to a rhombic crystal is continuous. 


C. Chiral symmetry breaking 

Finally, we investigate whether or not systems of equi¬ 
lateral and right-angled triangles show chiral symmetry 
breaking similar as was reported in the experiments of 
Ref. |9]. To this end, we perform Monte Carlo simulations 
of 5000 < N < 12000 triangles in the canonical ensem¬ 
ble. We calculate the orientational distribution function 
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Figure 7. Typical configurations of equilateral triangles at 
a packing fraction 0 = 0.97 (a) and 4> = 0.98 (b). The color 
coding of the particles is the same as in Fig. [6] Left-handed 
enantiomers are colored blue while right-handed enantiomers 
are colored red. The remaining particles are colored green. A 
clear phase boundary can be seen separating the two coexist¬ 
ing right- and left-handed chiral phases. 
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Figure 8. Close-up of the chiral phases at a packing fraction 
0 = 0.97. Top panel shows typical left-handed and right- 
handed chiral phases of equilateral triangles and the bot¬ 
tom panel displays the same for right-angled triangles. Left- 
handed enantiomers are colored blue while the right-handed 
enantiomers are colored red. 


P(0), where 6 is the angle that a triangle has with re¬ 
spect to a fixed axis (x-axis) as shown in Figs. 0b, c). 
Since the probability distribution to find anti-clockwise 
or clockwise orientational displacements should be sym¬ 
metric, i.e., P(6) = P(—0 ), we average the distributions 
for negative and positive 0 to get smoother probability 
distributions. We plot P(0) as a function of 0 in Fig.^d) 
for equilateral triangles and varying packing fractions <fi. 
We clearly observe that the unimodal distribution at low 
(j) splits into three distinct peaks at (j) x = 0.89 for equi¬ 
lateral triangles. The central peak corresponds to parti¬ 
cles oriented along the lattice vector while the remaining 
two peaks correspond to particles which have either anti¬ 
clockwise or clockwise orientational displacements. In 
Fig. 0 we show typical configurations for a system of 
equilateral triangles at = 0.97 and = 0.98. The tri- 



Figure 9. Projection of the center-of-masses of the equilateral 
triangles taken from 20 different equilibrium configurations at 
4> = 0.91. The hexagons are drawn to guide the eye to see 
the inherent honeycomb lattice. We have carved this picture 
from a larger system of N — 5000 for visual clarity. 

angles with a negative 6 , which are shifted anti-clockwise, 
are colored blue, whereas the triangles with a positive 6 
are colored red. The particles with an orientational dis¬ 
placement corresponding to the central peak in P(6) are 
colored green. Surprisingly, we find a clear phase sepa¬ 
ration between a phase with (blue) triangles that are ro¬ 
tated anti-clockwise and a phase with (red) particles that 
are twisted clockwise. The two coexisting phases are sep¬ 
arated by an interface of (green) particles that show no 
appreciable twist. We thus find an achiral triatic phase 
at 0 < 0.89, whereas the system phase separates into left- 
and right-handed chiral phases for > 0.89. Moreover, 
we find that the peaks corresponding to the two coex¬ 
isting chiral phases become more pronounced upon in¬ 
creasing 0, and hence the interfacial free energy increases 
with (j). We thus find that the phase behavior of hard tri¬ 
angles is remarkably similar as that of the Ising model, 
that shows at sufficiently low temperatures spontaneous 
magnetization and phase coexistence between two mag¬ 
netic phases. We therefore compared the orientational 
distribution function of the triangles with the probabil¬ 
ity distribution function of the magnetization of the Ising 
model in order to investigate if the demixing transition 
of triangles corresponds to the Ising universality class. In 
Fig. [lOj we show that the order parameter distribution 
functions do not match, and we conclude that the demix¬ 
ing transition of the enantiomers should correspond to 
another universality class, e.g., the six-state clock model. 
Finally, we wish to remark that the value of the most 
likely rotational shift 0 decreases upon increasing 0 as ex¬ 
pected since the rotational displacement equals zero for 
all triangles in the achiral crystal phase at close-packing. 

A similar chiral symmetry breaking and phase sep¬ 
aration is also observed for right-angled triangles (not 
shown). In this case, the transition from an achiral to a 
chiral phase occurs at <p x = 0.87. 

In Fig. [§J we show a close-up look of these chiral con¬ 
figurations for both the equilateral and right-angled tri¬ 
angles. We observe that the collective orientational dis- 
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Figure 10. The critical order parameter distribution function 
of the 2D Ising model P{m) as a function of the magneti¬ 
zation (black) and the orientational distribution function of 
triangles P{9) at a packing fraction cj) m 0.94. Both distri¬ 
bution functions are normalized to unit norm, variance and 
zero mean. The height of the orientational distribution of the 
triangles P(6) are rescaled to match P{m) 


placements of the triangles lead to a hexagonal lattice 
of clockwise or anti-clockwise chiral holes, which are sur¬ 
rounded by six triangles in the case of equilateral trian¬ 
gles. The appearance of these chiral holes due to the col¬ 
lective rotation of six triangles is also illustrated schemat¬ 
ically in Fig. (6^a). In the case of right-angled triangles 
the collective orientational displacements lead to a square 
lattice of chiral holes, which are surrounded by either 
four of eight triangles. We used curved arrows to in¬ 
dicate clockwise and anti-clockwise holes in Fig. [8j It 
is worth mentioning that for long simulation times the 
system should display either a pure left-handed or right- 
handed chiral phase, as it costs interfacial free energy to 
have a phase-separated configuration with an interface. 
Due to the long equilibration times the phase separated 
system frequently remains within the simulation times of 
our Monte Carlo runs. In order to investigate whether 
or not the ” twisted” triangles are still positioned on a 
regular lattice, we projected the center-of-masses of the 
equilateral triangles as obtained from 20 different equi¬ 
librated configurations at a packing fraction = 0.91 
on a plane in Fig. [9] We find that the center-of-masses 
of the particles form a regular honeycomb lattice with 
long-range positional order. 

In addition, we also computed the lateral shifts be¬ 
tween neighboring triangles at high densities as also com¬ 
puted by the authors of Refs. mm- Our results are in 
agreement with the earlier simulation results of Ref. HD 
that there is no split in the probability distributions of 
these lateral shifts for perfect hard triangles. 


II. PHASE DIAGRAM AND CONCLUSIONS 

In summary, we have studied a two-dimensional system 
of equilateral triangles and right-angled isosceles trian¬ 
gles using large-scale Monte Carlo simulations. We have 
computed the equations of state, and bond-orientational 
and molecular orientational order parameters as a func¬ 
tion of packing fraction 0. In addition, we calculated 
the free energies as a function of packing fraction for the 
isotropic fluid phase, the liquid crystal phase, and solid 
phase. We also measured the spatial correlation functions 
for the translational, bond-orientational, and molecular 
orientational order. We mapped out the phase diagram 
of both equilateral triangles and right-angled triangles by 
combining these results. In Figs. |3je,f) we summarized 
the phase behavior using different colors. We indicate 
the different phase transitions by vertical dotted lines as 
a guide to the eye across the different graphs. We show 
that hard equilateral triangles and hard right-angled tri¬ 
angles undergo a phase transition from an isotropic phase 
to a triatic and rhombic liquid crystal phase, respec¬ 
tively. The phase transition from the isotropic to triatic 
liquid crystal phase is continuous for equilateral trian¬ 
gles, whereas we find a first-order phase transition from 
the isotropic fluid to the rhombic liquid crystal phase 
for the right-angled triangles with a coexistence region 
4> £ [0.733,0.782]. With increasing pressure these liquid 
crystal phases continuously transform to their respective 
close-packed crystal structures. These close-packed crys¬ 
talline phases exhibit at sufficiently high packing frac¬ 
tions spontaneous chiral symmetry breaking as the trian¬ 
gles rotate either in clockwise or anti-clockwise direction 
with respect to a fixed lattice vector. We denote the chi¬ 
ral triatic phase and the chiral rhombic phase by T x and 
R x , respectively, in the phase diagram of Figs. [3](e,f ). We 
also observe a spontaneous purely entropy-driven demix¬ 
ing of the ” enantiomers” resulting in phase coexistence 
of the left- and right-handed chiral phase with a clear in¬ 
terface. To the best of our knowledge, our work presents 
the first observation of a spontaneous macroscopic chi¬ 
ral symmetry breaking and entropy-driven demixing of 
” enantiomers” of achiral building blocks. Finally, we 
wish to remark that the isotropic-to-liquid-crystal phase 
transition point in equilateral triangles as determined in 
experiments and in an earlier simulation study mm are 
15% off from our simulation results. Additionally, the 
EOS as shown in Fig. 1 of Ref. [12] does not match with 
our EOS obtained from our isotensic NPT Monte Carlo 
simulations. We attribute this discrepancy with earlier 
simulation results m to the fact that these molecular 
dynamics simulations were performed with a fixed box 
shape, which may lead to non-zero stress. We verified 
this by Monte Carlo simulations of hard triangles in a 
fixed box shape, which indeed show that the isotropic- 
to-liquid-crystal phase transition happens at lower pack¬ 
ing fraction compared to simulations with a variable box 
shape. The mismatch with the experimental |9j isotropic- 
liquid-crystal phase transition point is likely due to the 
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fact that the particle interactions in the experimental 
system cannot be described by excluded-volume interac¬ 


tions, which may be caused by the presence of depletants, 
charges, and polydispersity |9l [23]. 
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